Quantum kinetic theory VII: The influence of vapor dynamics on condensate growth 
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We extend earlier models of the growth of a Bose-Einstein condensate [|l|-^|] to include the full dynamical 
effects of the thermal cloud by numerically solving a modified quantum Boltzmann equation. We determine 
the regime in which the assumptions of the simple model of [^] are a reasonable approximation, and compare 
our new results with those that were earlier compared with experimental data. We find good agreement with 
our earlier modelling, except at higher condensate fractions, for which a significant speedup is found. We also 
investigate the effect of the final temperature on condensate growth, and find that this has a surprisingly small 
effect. 

The particular discrepancy between theory and experiment found in our earlier model remains, since the 
speedup found in these computations does not occur in the parameter regime specified in the experiment. 



I. INTRODUCTION 

The fundamental process in the growth of a Bose-Einstein 
condensate is that of bosonic stimulation, by which atoms are 
scattered into and out of the condensate at rates enhanced by 
a factor proportional to the number of atoms in the conden- 
sate. This was first quantitatively considered by Gardiner et al. 
[jjj], in a paper which treated the idealized case of the growth 
of a condensate from an nondepletable "bath" of atoms at a 
fixed positive chemical potential /i and temperature T. This 
gave rise to a simple and elegant formula known as the simple 
growth equation 

h = 2W+(n ) Ml- e [Mc(no)-M]/k^ Uq + A ; (1) 

in which no is the population of the condensate, /i is the chem- 
ical potential of the thermal cloud, and /Uc(^o) is the con- 
densate eigenvalue. The prefactor W + (no) is a rate with an 
expression derived from quantum kinetic theory [0], which 
was estimated approximately in ^ by using a classical Boltz- 
mann distribution. To go beyond the Boltzmann approxima- 
tion for W + involves a very much more detailed treatment 
of the populations of the trap levels with energy less than fx, 
since the equilibrium Bose-Einstein distribution for fx > is 
not consistent with energies less than fi. In other words, the 
populations of the lower trap levels cannot be treated as time- 
independent, and thus the dynamics of growth must include at 
least this range of trap levels as well as the condensate level. 
Therefore in [||3|] we considered a less simplified model, cov- 
ering a range of energies up to a cut-off, Er, above which the 
system was assumed to be a thermal cloud with a fixed tem- 
perature and chemical potential. Equations were derived for 
the rate of growth of these levels along with the condensate, 
and the rates at which particles from the thermal bath scat- 
tered these quasi-particles between levels within the conden- 



sate band. The results of calculations showed that a speedup 
of the growth rate by a factor of the order of 3^1 compared 
to the simple growth equation could be expected, and that the 
initial part of the growth curve would be modified, leading to a 
much sharper onset of the initiation of the condensate growth. 

The only experiment that has been done on condensate 
growth [[FT]] was then under way. In these experiments clouds 
of sodium atoms were cooled to just above the transition tem- 
perature, at which point the high energy tail of the distribution 
was rapidly removed by an very severe RF "cut", where the 
frequency of the RF field was quickly ramped down. After a 
short period of equilibration, the resulting vapor distribution 
was found to be similar to the assumptions of our theoretical 
treatments, and condensate growth followed promptly. The 
results obtained were fitted to solutions of the simple growth 
equation ([j]). When experimental results became available a 
speedup of about the predicted factor was found, and indeed 
the higher temperature results agreed very well with the theo- 
retical predictions. At lower temperatures there was still some 
disparity; the theory predicted a slower rate of growth with 
decreasing temperature, but experimentally the opposite was 
observed. 

The situation in which we now find ourselves leaves no al- 
ternative other than to address the remaining approximations. 
In our previous work we have made four major approxima- 
tions: 

(i) The part of the vapor with energies higher than Er has 
been treated as being time independent. 

(ii) The energy levels above the condensate level were mod- 
ified phenomenologically to account for the fact that 
they must always be greater than the condensate chem- 
ical potential, which rises as the condensate grows. 

(iii) We treated all levels as being particle-like, on the 
grounds that detailed calculations [|j~2j] have shown that 



1 



only a very small proportion of excitations of a trapped 
Bose gas are not of this kind. 

(iv) We have used the quantum Boltzmann equation in a er- 
godic form, in which all levels of a similar energy are 
assumed to be equally occupied. 

In this paper we will no longer require the first two of these 
approximations. Abandoning the first means that we are re- 
quired to take care of all kinds of collisions which can occur, 
and thus treat the time-dependence of all levels. This comes 
at a dramatic increase in both the computation time required 
(hours rather than seconds) and the precision of algorithms 
required. We also use a density of states which should be 
close to the actual density of states as the condensate grows, 
thereby avoiding the phenomenological modification of en- 
ergy levels. However, we still treat all of the levels as being 
particle-like, since it seems unlikely that the few non-particle- 
like excitations will have a significant effect on the growth 
as a whole. The ergodic form of the quantum Boltzmann is 
needed to make the computations tractable, and is of neces- 
sity retained. 



H. FORMALISM 

The basis of our method is Quantum Kinetic theory, a full 
exposition of which is given in Ref. M]. These develop a com- 
plete framework for the study of a trapped Bose gas in a set of 
master equations. The full solution of these equations is not 
feasible, however, and therefore some type of approximation 
must be made. The basic structure of the method used here 
is essentially the same as that of QKVI, the major difference 
being that all time-dependence of the distribution function is 
retained. As explained in ^[u}, quantum kinetic theory leads 
to a model which can be viewed as a modification of the quan- 
tum Boltzmann equation in which 

i) The condensate wavefunction and energy eigenvalue — 
the condensate chemical potential /xc(^o) — are given 
by the solution of the time-independent Gross- 
Pitaevskii equation with no atoms. 

ii) The trap levels above the condensate level are the quasi- 
particle levels appropriate to the condensate wavefunc- 
tion. This leads to a density of states for the trap levels 
whic h is s ubstantially modified, as discussed below in 
Sect. Mb 



iii) The transfer of atoms between levels is given by a mod- 
ified quantum Boltzmann equation (MQBE) in the en- 
ergy representation. This makes the ergodic assump- 
tion; that the distribution function depends only on en- 
ergy. 



A. The ergodic form of the quantum Boltzmann equation 

The derivation of the ergodic form of the quantum Boltz- 
mann equation used by [|l3]| is particular to the undeformed 



harmonic potential, and we give here a derivation appropriate 
to our case, in which the density of states can change with 
time as the condensate grows. We bin the phase space into 
energy bands labeled by the index n with energies in a range 

D n {t) = (e n ( t )-^,e n (t) + of width 5s n (t), 

and these widths change in time so that the number of states 
within each bin, g n is constant in time. 

Starting from the full quantum Boltzmann equation the er- 
godic approximation is expressed in terms of this binned de- 
scription as follows: We set /(x, K, t) equal to a constant, 
f n , when e(x, K, t) = h 2 K. 2 /2m + V c ff(x,t) is inside the 
nth bin, i.e., e(x, K, t) g D n (t). (Here T4fj(x, t) is the poten- 
tial of the trap, as modified by the mean field arising from the 
presence of the condensate, as explained in QKV and QKVI.) 

Thus we can approximate 



a/(x,K,t) 

dt 



% if e(x,K,t)eD n (t), (2) 
at 



In order to derive the ergodic quantum Boltzmann equation, 
we define the indicator function x„(x, K, i) of the nth bin 

£>„(<) by 

X„(x,K,i) = l if e(x,K,t)eD n (t), (3) 
= otherwise. (4) 

The number of states in the bin n will be given by g n = 
J d 3 xd 3 K x„(x, K, t) /h 3 , and is held fixed. 

The formal statement of the binned approximation is 



/(x,K,i)^^./-„ Xll (x,K,i), 



(5) 



and the ergodic quantum Boltzmann equation is derived by 
substituting (Q) into the various parts of the quantum Boltz- 
mann equation as follows. For the time derivative part we 
make this replacement, and project onto D n (t), getting 



d 3 xd 3 K df(x,K,t) 
X„(x,K,i)- 



dt 



din 

dt 



9n-ET- ( 6 ) 



[Note that the expansion (Q) would mean that delta func- 
tion singularities at the upper and lower boundaries of D„ (t) 
would arise by differentiating /(x, K,t) as defined in (g), 
but the condition that g n be fixed means that these are of 
equal and opposite weight, and cancel when integrated over 
D n (t), givng a result consistent with (p)).] We now replace 
9/(x, K, t) /dt on the left hand side of (^ by the collison inte- 
gral that appears on the right hand side of the quantum Boltz- 
mann equation, and substitute for /(x, K,t) in the collision 
integral using (ph. [The streaming terms give no contribution, 
since the form (ph is a function of the energy e(x, K, t).] 

This leads to the ergodic quantum Boltzmann equation in 
the form 
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9n^ = ^ E ifpM 1 + /r)(l + /n) - (1 + /p)(1 + /,)/r/«} / / f d 3 K 2 f d 3 K 3 

771 pqr •* J J J 

x x „(x, K!,t)x P (x, K 2 , t) x ,(x, K 3 , t)xr(x, K, «)* (K x + K 2 - K 3 - K) 

x5 (e(x, Ki, t) + e(x, K 2 , t) - e(x, K 3 , t) - e(x, K, *)) . (7) 
The final integral is now approximated by the method of Jl3| ] to give [a and ui are defined in Sect. Ill | 

9nj£ = E {/p/^ 1 + /r)(l + /n) - (1 + U)0- + fg)frfn} M(p, g, r, »)A(p, g, r, n). (8) 

pqr 

Here A(p, g, r, n) is a function which expresses the overall energy conservation, and is defined by 

a / \ -. i I . \Se p + 6e q + 6s r + 6s n \ 

A(p,q,r,n) = 1 when \e p + e q -e r - e n \ < — ^— , (9) 

= otherwise. (10) 



Because we approximate /(x, K, t) by a constant value 
within each D n (t), energy conservation means that E = 
Yin £ ngnfn(t) is constant. This follows from energy conser- 
vation in the full quantum Boltzmann equation, which also 
implies that 

E A (P' r ' n ) M (P> 1, r, n)(e r + e„) 

rn 

= {e p + e q )^2A(p,q,r,n)M(p,q,r,n) (11) 

This is the limit to which the binning procedure defines energy 
conservation. 
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FIG. 1. Qualitative picture of the compression of the quantum 
levels above the condensate mode as the condensate eigenvalue in- 
creases. 



results in a compression in energy space of the quantum lev- 
els immediately above the condensate, (see Fig. |l|) and has an 
important effect on the evolution of the cloud. 

The correct description of the quantum levels immediately 
above the ground state when there is a significant conden- 
sate population requires a quasiparticle transformation. This 
is computationally difficult, however, so we make use of a 
single-particle approximation for these states. This should 
be reasonable, as most of the growth dynamics will involve 
higher lying states that will be almost unaffected by the pres- 
ence of the condensate. In [j^] we did this using a linear inter- 
polation of the density of states; here we use an approximate 
treatment based on the Thomas-Fermi approximation. 



A. Condensate chemical potential /ic(«o) 

We consider a harmonic trap with a geometric mean fre- 
quency of 

U = (uj x uj y uj z ) 1/3 . (12) 

We include the mean field effects via a Thomas-Fermi approx- 
imation for the condensate eigenvalue, which is directly re- 
lated to the number of atoms in the condensate mode. As in 
[^.[[(JL we use a modified form of this relation in order to give 
a smooth transition to the correct harmonic oscillator value 
when the condensate number is small: 



III. DETAILS OF MODEL 

The most important aspect of our model is the inclusion of 
the mean field effects of the condensate. As the population 
of the condensate increases, the absolute energy of the con- 
densate level also rises due to the atomic interactions. This 



McC^o) = ct no + (3hLu/2a) 5 ' 2 



2/5 



(13) 



where a = (15aa;m 1/ ' 2 ?i 2 /4v / 2) 2 ^ 5 , and a is the atomic s- 
wave scattering length. Thus, for no = we have nc (0) = 
£ = 3hu>/2. 
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B. Density of states g(e) 



Q], which leads to a density of states of the form 



We assume a single particle energy spectrum with a 
Bogoliubov-like dispersion relation, as in Timmermans et al. 



4 Mc(np) s 

7T (huj) 3 



; r- — 1 I / d,X\/l — X 

Mc(n ) / Jo 



#/wS)-i) 2 + i 2 - 



This integral can be carried out analytically; the result is 

-2 f 



g(s,n ) 



I - /- Mc(n ) 



2{T%ujf 



1 + (/*c(no)/e)+ 1- 



</2 



1/2 



V / (e/Atc(n )-l) 2 + a: 2 



e/mc(".o) 



dxy/x 



Mc(n ) 



(14) 



£/Vc(™o) - 1, 



(15) 



where 



<7i(x) = - [VaVl - s(l - 2x) - sin 1 (v / ^)] , 



92 (x) 



4V2 



/— , /l+X + V2^\ f 7T . _! / 1-1 

'2x + x log I , ) - <! - + sin 1 



VT+x^ 



(16) 
(17) 



This is plotted in Fig. |2|, along with that for the ideal gas. Thus 
the density of states of the trap varies smoothly as the conden- 
sate grows. 




FIG. 2. The modified density of states (solid curve) compared 
with non-interacting function (dashed curve) for the harmonic trap. 



IV. NUMERICAL METHODS 



A. Representation 



The bins we shall choose for the representation of the dis- 
tribution in terms of the quantities /„ as in (|5|) are divided into 
two distinct regions, as shown diagrammatically in Fig. ||. The 



lowest energy region corresponds essentially to the conden- 
sate band Rc of [ffM^,|o| . This is the region in which /„ is 
rapidly varying in the regime of quantum degeneracy, and is 
described by a series of fine-grained energy bins up to an en- 
ergy En sa 3^0(^0, max)- The condensate is a single quantum 
state represented by the lowest energy bin. 
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FIG. 3. The numerical representation of the system with a con- 
densate of 2.3 x 10 6 atoms at a temperature of 590nK. Rc is 
the condensate band, which is fine-grained, whereas Rnc is the 
non-condensate band, which is coarse-grained. The division between 
the two bands is fixed at Er. The condensate energy is derived from 
the Thomas-Fermi approximation. 

As the number of particles in the condensate changes, 
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the energy of the condensate level c hang es according to the 
Thomas-Fermi approximation of Eq. [lj. Thus the total en- 
ergy width of Rc decreases as the condensate grows. 

We represent Rc by a fixed number of energy bins of equal 
width 5e n with a midpoint of e n . 

As the condensate energy increases, we adjust e n and fe ra 
between integration timesteps, such that the widths of the bins 
remain equal. This is done by redistributing the numbers of 
particles into new bins after each timestep, and thus does not 
contradict the requirement that g n is fixed during the timestep. 
We find that this is the most simple procedure for the calcula- 
tion of rates in and out of these levels. We choose the number 
of bins to be sufficient such that the width is not more than 
about Se n ~ 5fuD. 

The high energy region corresponds to the thermal bath 
of our previous papers. This is the region in which /„ is 
slowly varying, and therefore the energy bins are considerably 
broader (up to QATiu in the results presented in this paper). 
The evaporative cooling is carried out by the sudden removal 
of population of the bins in this region with e„ > e cut . 



B. Solution 

There are four different types of collision that can occur 
given our numerical description of the system. These are de- 
picted in Fig. ^. 




(a) Growth 




(b) Scattering 



condensate band (along with the reverse process). 

(b) Scattering: A particle in Rnc collides with a particle 
in the condensate band, with one particle remaining in 
Rc- 

(c) Internal: Two particles within the condensate band col- 
lide with at least one of these particles remaining in 
Rnc a f ter the collision. 

(d) Thermal: This involves all particles involved in the col- 
lision coming from the non-condensate band and re- 
maining there. 

Our first description of condensate growth [[[]] considered only 
process (a). The next calculation involved both process 
(a) and (b). The calculations presented below include all four 
processes, allowing us to determine whether the earlier ap- 
proximations were justified. 

The computation of the rates of processes (a) and (b) is 
made difficult because of the different energy scales of the two 
regions of the distribution function. Our solution is to inter- 
polate the distribution function /„ in Rmc (non-condensate 
band) such that the bin sizes are reduced to be the same as for 
Rc (the condensate band). The rates are then calculated using 
this interpolated distribution function, now consisting of more 
than one thousand bins, and the rates for the large bins of the 
non-condensate band are found by summing the rates of the 
appropriate interpolated bins. 

We have found that these rates are extremely sensitive to the 
accuracy of the numerical interpolation — small errors lead to 
inconsistencies in the solutions of the MQBE. This procedure 
is more efficient than simply using the same bin size for the 
whole distribution, as there are only a small number of bins 
for the condensate band. 



C. Algorithm 

The algorithm we use to solve the MQBE is summarised as 
follows 




(c) Internal 




(d) Thermal 



FIG. 4. The four different collision types that can occur in our 
numerical description. 

(a) Growth: This involves two particles in Rnc colliding, 
resulting in the transfer of one of the particles to the 



(1) Calculate the collision summation for all types of col- 
lisions, keeping the density of states, and the energies 
of the levels in the condensate band Rc fixed. The 
distribution function f n (t + 5t) is calculated using an 
embedded 4th order Runge-Kutta method, using Cash- 
Karp parameters [ |T5| ] . 

(2) The quantity M(p, g, r, n) defined by (^|) expresses all 
the overlap integrals, and is quite difficult to compute 
exactly. In our computations we have simply set this to 
correspond to the value found in [pj|], i.e., we set 



M(p, q, r, n) = g m i a ( P ,q,r,n), 



(18) 



and express energy conservation in a simplified form, 
using the fact that the energy bins will be chosen equally 
spaced, by choosing a Kronecker delta form 
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A(p, q,r,n) — > 8{p + q,r + n). 



(19) 



The difference between these two forms clearly goes to 
zero as the bins become very narrow. 

It has been explicitly checked that in practice energy is 
conserved to a very high degree of accuracy throughout 
the calculation. 

(3) As a result of the time step, the condensate population 
will have changed. This causes the density of states to 
alter slightly, along with the positions and widths of the 
energy bins in the condensate band, as all these quanti- 
ties are determine d by the condensate number no. The 
derivation in Sect. II A shows that the populations g n f n 
are of the bins which move with the change of energy 
levels and density of states as the condensate grows so 
as to maintain the number of levels g n in the bin con- 
stant. Therefore after the Runge-Kutta timestep, the 
numbers g„fn represent the numbers of particles in the 
bins determined by the appropriate energy levels after 
that step. 

(4) As a result of the preceding step, the bins will no longer 
of equal width, so we rebin the numbers of atoms into a 
new set of equally spaced bins, as explained Sect. IV A 



To ensure total number conservation of particles, we 
keep the number of particles in each bin, g„f n constant 
when we adjust the energies and widths of the bins. As 
the change in the density of states and the width of each 
bin is determined by the condensate number, the oc- 
cupation per energy level of the nth bin, f n , must be 
altered slightly to ensure number conservation. 

(5) We now continue with step (1). 

The change in nc(no) with each time step, and hence the 
shifts in the energy of the bins in Rc is very small. Hence, the 
adjustment of the distribution function due to step (3) is tiny, 
much smaller than the change due to step (1). 

The method has been tested by altering the position of Er 
and width of the energy bins of Rnc an d Rc- We have found 
that the solution is independent of the value of Er over a large 
range of values of these parameters. 



V. RESULTS 

In this paper we present the results of simulations mod- 
elling the experiments described in [ [TT| ] . In these experiments 
a cloud of sodium atoms confined in a 'cigar' shaped magnetic 
trap was evaporatively cooled to just above the Bose-Einstein 
transition temperature. Then, in a period of 10ms the high en- 
ergy tail of the distribution was removed with a very rapid and 
rather severe RF cut. The condensate was then manifested by 
the formation of a sharp peak in the density distribution. 

We have carried out a full investigation of the effect of vary- 
ing the initial cloud parameters has on the growth of the con- 
densate for the trap configuration described in [pTJ]. In this 



paper we concentrate on a comparison of these results with 
our earlier theoretical model. To model these experiments, we 
begin our simulations with an equilibrium Bose-Einstein dis- 
tribution, with temperature and chemical potential fj,- m i t and 
truncate it at an energy e cut = rjkTi, which represents the sys- 
tem at the end of the RF sweep. This is then allowed to evolve 
in time, until the gas once again approaches an equilibrium 
the appropriate equilibrium Bose-Einstein in the presence of 
a condensate. This is pictured schematically in Fig. g. 
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FIG. 5. Snapshots of the distribution function for a simulation 
with initial conditions /iinit = — WOhui, Tj = 1119nK, and 
n = 2.83. This results in a condensate with no = 7.5 x 10 6 atoms 
at a temperature of T/ = 830nK. For clarity, the condensate itself is 
not depicted, but the presence of a significant amount of condensate 
has the effect of displacing the left hand ends of the curves (d)-(f) 
an amount /ic(no)/fcT; from the axis. The growth curve for this 
simulation is shown in Fig. |^(a). 



Because of the ergodic assumption, the MQBE that we sim- 
ulate depends only on the geometric average of the trapping 
frequencies ui = (u) x ujyUJ z ) 1 /' i . There is likely to be some 
type of experimental dependence on the actual trap geome- 
try which is not included in our simulation; however in the 
regime kT ^> %ui this should be small. The trap parameters 
of Q were (u x , u y , w z ) = 2tt x (82.3, 82.3, 18)Hz , giving 
Cu = 2tt x 49.6Hz. 
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A. Matching the experimental data 



The main source of quantitative experimental data of con- 
densate growth generally available is Fig. 5 of [|j"l|]. This gives 
growth rates as a function of final condensate number and tem- 
perature rather than the initial conditions. Whereas the growth 
curves calculated in [|lj|^] required these parameters as inputs, 
the calculations presented here require three different input 
parameters; the initial number of atoms in the system Ni (and 
hence the initial chemical potential /imitX the initial tempera- 
ture Tj, and the position of the cut energy rjfcTj. 



the measured final temperature Tf, we can calculate the to- 
tal energy E tot and number of atoms iV tot in the system at 
the end of the experiment, completely characterizing the final 
state of the gas. 



iVtot = n Q + 

e„>)ic(no) 
oo 

E tot =E (n a )+ 



On 



exp[(e„ - nc{n Q ))/kTf\ - 1 ' 



£n>McOo) 



(20) 



exp[(e„ - nc{m))/kT f ] - 1 ' 



(21) 
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FIG. 6. Growth of a condensate with no = 7.5 x 10 6 , 
Tf — 830nK. Solid lines /Uinit = 0, dotted lines nuxn = — 40to, 
dashed lines fii n it = — 100 to. (a) Population of condensate versus 
time. Grey curve is the solution for model of |Q|, (b) Chemical po- 
tential nc(no) of condensate (lower curves) and effective chemical 
potential (i c ff of thermal cloud (upper curves). 



Given the final parameters supplied in [[UJ], it is possible 
to calculate a set of initial conditions that we require. As we 
know the final condensate number, we can calculate the value 
of the chemical potential of the gas using the Thomas-Fermi 
approximation for the condensate eigenvalue, Eq. |l3|. This 
gives a density of states according to Eq. [l5|, and along with 



We now want to find an initial distribution that would have 
the same total energy and number of atoms if truncated at 
£cut = T]kTi. If we specify an initial chemical potential for 
the distribution /ii n it, we can self-consistently solve for the 
parameters Ti and i] from the following non-linear set of equa- 
tions 



JjfcT; 

N tot = £ 

e„=3fiu)/2 
rikTi 



9n 

exp[(e„ - /Uinit) /kTi] - 1 : 



(22) 



(23) 



This gives the input parameters for our simulation, and we 
can now calculate growth curves starting with initially differ- 
ent clouds, but resulting in the same final condensate number 
and temperature. 



B. Typical results 

A sample set of growth curves is presented in Fig. |^(a), for 
a condensate with 7.5 x 10 6 atoms at a final temperature of 
830nK and a condensate fraction of 10.4%. The initial pa- 
rameters for the curves are given in Table [I| 

As can be seen the curves are very similar, and arguably 
would be difficult to distinguish in experiment. The main 
difference is the further the system starts from the transition 
point (i.e. the more negative the initial chemical potential), 
the longer the initiation time but the steeper the growth curve. 
This trend continues as /ij n i t becomes more negative. 



/iinit(to) 


T (nK) 


A^IO 6 ) 


V 


£ cut (huj) 





1000 


89.1 


3.82 


1605 


-40 


1080 


100.1 


3.31 


1503 


-100 


1119 


117.6 


2.83 


1419 



TABLE I. Parameters for the formation of a condensate with 



7.5 x 10 6 atoms at a temperature of Tf = 830nK from an 
uncondensed thermal cloud. The growth curves are plotted in Fig. ^. 
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1. Effective chemical potential 
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To facilitate understanding of these results, we introduce 
the concept of an effective chemical potential /i e g for the non- 
condensate band. We do this by fitting a Bose-Einstein distri- 
bution to the lowest energy bins of Rnc as a function of time. 
Obviously, the chemical potential is undefined when the sys- 
tem is not in equilibrium, but as has been noted for the clas- 
sical Boltzmann equation Qlq], the distribution function tends 
to resemble an equilibrium distribution as evaporative cooling 
proceeds. The effective chemical potential is not unique — it is 
dependent on the particular choice of the energy cutoff Er. It 
gives a good indication of the "state" of the non-condensate, 
however, since the majority of the particles entering the con- 
densate after a collision come from these bins. In this paper 
/j, e ff was computed by a linear fit to log[l + 1/ /„] of the first 
ten bins of the noncondensate band, with the intercept giving 
/j, e g and the gradient the temperature. 



2. Interpretation 



We find that all the results presented in this paper can be 
qualitatively understood in terms of the simple growth equa- 
tion Eq. with the vapour chemical potential fx replaced by 
the effective chemical potential /i e g of the thermal cloud. 

The simple growth equation requires fj, e g > /ic(^o) f° r 
condensate growth to occur. In Fig. ||(b) we plot the effective 
chemical potential fi c g of the thermal cloud and the chem- 
ical potential of condensate Hc(no)- This graphs helps ex- 
plain the two effects noted above — longer initiation time and 
a steeper growth curve for the fj,^ = — 100?wl> case. Firstly, 
the inversion of the chemical potentials for this simulation oc- 
curs at a later time than for fi ull t = 0, causing the stimu- 
lated growth to begin later. This is because the initial cloud 
for the /iinit = — lOOto simulation is further from the tran- 
sition point at t = 0. Secondly, the effective chemical po- 
tential of the thermal cloud rises more steeply, meaning that 
A^cff — Mc(%) is larger, and therefore the rate of condensate 
growth is increased. 



C. Comparison with earlier model 



In Fig. ||(a) we have also plotted the growth curve that is 
calculated for these final condensate parameters by the model 
of [[3J, which we refer to as the simple model (not to be mis- 
taken with the solution of the simple growth equation [l]). For 
this earlier model the initial condensate number is indetermi- 
nate, whereas for the detailed calculation presented here the 
initial distribution is Bose-Einstein, with the zero of the time 
axis being the removal of the high energy tail. 
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FIG. 7 

sate fraction of 24.1%, no = 7.5 
/iinit = 0, dashed lines fj,- in it = 



densate versus time. Grey line is the solution for model of Pfl. (b) 
Chemical potential of condensate (lower curves) and effective chem- 
ical potential of thermal cloud (upper curves). 

For these particular parameters, it turns out that the results 
of the full calculation of the growth curve give very simi- 
lar results to the previous model, with the initial condensate 
number adjusted appropriately. This is not surprising; indeed, 
from Fig. ||(b) we can see that the approximation of the ther- 
mal cloud by a constant chemical potential (i.e. the cloud is 
not depleted) is good for the region where the condensate be- 
comes macroscopic. 

For larger condensate fractions, however, the principal con- 
dition assumed in the model of [j^, that the chemical potential 
of the vapor can be treated as approximately constant, is no 
longer satisfied. In Fig. [7(a) we plot the growth of the same 
size condensate as in Fig. 5[ (that is, 7.5 x 10 6 atoms) but at a 
lower final temperature of 590nK. In this situation the conden- 
sate fraction increases to 24. 1 %, and so there is considerable 
depletion of the thermal cloud. The effect of this can be seen 
in Fig. 0(b). The difference between the vapor and conden- 
sate chemical potentials /i e ff — /Uc(fto) initially increases to 
much larger values than for the simple model, where /x e ff is 
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held constant at its final equilibrium value. It is this fact that 
causes more rapid growth. 

As the condensate continues to grow, it begins to signif- 
icantly deplete the thermal cloud, causing /i e g to decrease 
from its maximum. It is the "overshoot" of /x e g from the fi- 
nal equilibrium value that the model of [[|] and QKVI cannot 
take account of. This overshoot only occurs for final conden- 
sate fractions of more than about ten percent; hence up to this 
value the simple model should be sufficient. 
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D. Effect of the final temperature on condensate growth 



We have investigated the effect that final temperature has on 
the growth of a condensate of a fixed number. All simulations 
were begun with /ij n i t = 0, since the initial chemical potential 
has little effect on the overall shape of the growth curves. This 
determines the other parameters, Tj and r) and the initial con- 
ditions are shown in Table [n| The results of these simulations 
are presented in Fig. [| 

We find the somewhat surprising result that the growth 
curves do not change significantly over a very large tempera- 
ture range for the same size condensate. In fact, a condensate 
formed at 600nK grows more slowly than at 400nK for these 
parameters. As the temperature is increased further, however, 
the growth rate increases again, and at a final temperature of 
1 fiK the growth rate is faster than at 400nK. This effect has 
also been observed for both larger (7.5 x 10 6 ) and smaller 
(1 x 10 6 ) condensates. 

This observation can once again be interpreted using the 
simple growth equation Although W + (no) increases with 
temperature (approximately as T 2 as shown in [[[J), the maxi- 
mum value of u c ff — /ic( n o) achieved via evaporative cooling 
decreases with temperature for a fixed condensate number, as 
the cut required is less severe and the final condensate frac- 
tion is smaller. Also, the term in the curly brackets of Eq. [j] 
is approximately proportional to T -1 for most regimes. The 
end result is that the decrease in this term compensates for the 
increase in W + (no), giving growth curves that are very simi- 
lar for the different simulations. Once the "overshoot" of the 
thermal cloud chemical potential ceases to occur (when the 
evaporative cooling cut is not as severe), the growth rate be- 
gins to increase with temperature as predicted by the model of 
M\ and QKVI. 
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Time (s) 

FIG. 8. Growth of a condensate with a final condensate size of 
2.5 x 10 6 atoms from a vapor with /iinit = 0. The dotted line is 
for a final temperature of 400nK, dashed 600nK, and solid 1/iK. (a) 
Growth curves, (b) Chemical potential of condensate (lower curves) 
and thermal cloud (upper curves). 



T>(nK) 


Ti(nK) 


iV s (10 6 ) 


V 




Condensate fraction 


400 
600 
1000 


622.0 
707.3 
1064.8 


21.5 
31.6 
107.7 


2.19 
4.03 
5.87 


572 
1198 
2629 


0.253 
0.099 
0.025 



TABLE II. Parameters for the formation of condensate with 
no = 2.5 x 10 6 atoms from an uncondensed thermal cloud with 
Minit = 0. The growth curves are presented in Fig. ||. 
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E. Effect of size on condensate growth 
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FIG. 9. Growth of a condensate with at a final temperature of 
590nK, starting from an uncondensed thermal cloud with fi ln i t = 0. 
Solid line 7.5 x 10 6 atoms, dotted line 5.0 x 10 6 atoms, dashed 
2.3 x 10 6 atoms. The dashed line is for the same parameters as the 
lower temperature curves in (a) Growth curves, with the solu- 
tion to the model of |^] in grey, (b) Chemical potential of condensate 
(lower curves) and thermal cloud (upper curves). 



n (10 6 ) 


Ti(nK) 


A^IO 6 ) 


V 


£ cut (tlL0) 


Condensate fraction 


2.3 


692.5 


29.6 


4.07 


1186 


0.095 


5.0 


794.6 


44.7 


2.91 


973 


0.179 


7.5 


897.9 


64.6 


2.29 


865 


0.239 



TABLE III. Parameters for the formation of condensates at T/ 
= 590nK from an uncondensed thermal cloud with /tihiit = 0. The 
growth curves are presented in Fig. M 



Finally we have performed some simulations of the forma- 
tion of a condensate at a fixed final temperature, but of a vary- 
ing size. The parameters for these simulations are given in 
Table H and the growth curves are plotted in Fig. [ja). We 
find that the larger the condensate, the more rapidly it grows. 
The initial clouds required to form the larger condensates not 
only start at a higher temperature (and thus have a higher col- 
lision rate to begin with), but also they need to be truncated 
more severely, causing a larger difference in the chemical po- 
tentials, as seen in Fig. ||(b). Thus instead of these effects 
negating each other as in the previous section, here they tend 
to reinforce one another. This causes the growth rate to be 
highly sensitive to the final number of atoms in the conden- 
sate for a fixed final temperature. 

For further comparison with the previous model, in 
Fig. ^(a) the dashed curve is for the same parameters as for 
the lower temperature results of [gj], whose prediction is plot- 
ted in grey — as can be seen, the two methods are in very good 
agreement with each other for this choice of parameters. It is 
this particular set of final parameters for which the discrep- 
ancy between theory and experiment remains. 



F. The appropriate choice of parameters 

In our computations we have taken some care to make sure 
that we can give our results as a function of the experimen- 
tally measured final temperature Tf and condensate number 
no- Nevertheless, it can be seen from our results that this 
can give rise to counterintuitive behavior, such as the fact that 
under the condition of a given final condensate number, the 
growth rate seems to be largely independent of temperature, 

This effect 



because of the cancellation noted in Sect. VD. 
has its origin in the quite simple fact that with a sufficiently 
severe cut it is impossible to separate the process of equilibra- 
tion of the vapor distribution to a quasiequilibrium from the 
actual process of growth of the condensate. In other words, the 
attempt to implement the "ideal" experiment in which a con- 
densate grows from a vapor with a constant chemical potential 
and temperature cannot succeed with a sufficiently large cut. 
Under these conditions, the initial temperature differs quite 
strongly from the final temperature, and as well, the number 
of atoms required to produce the condensate is so large that 
the vapor cannot be characterized by a slowly varying chemi- 
cal potential during most of the growth process. 



G. Comparison with experiment 



1. Comparison with MIT fits 



The most quantitative data available from [ |1 1| ] is in their 
Fig. 5, in which results are presented as parameters extracted 
from fits to the simple growth equation |IJ In ^ we took 
two clusters of data from this figure, at the extremes of the 
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temperature range for which measurements were made, and 
compared the theoretical results with the fitted experimental 
curves. At the higher temperature of 830nK the results were 
in good agreement with experiment, but at 590nK they dif- 
fered significantly, the experimental growth rate being about 
three times faster than the theoretical result. 

We have performed the same calculations using the detailed 
model. The results for 830nK arepresented in Fig. ||and those 
for 590nK are presented in Fig. |. There is a good match be- 
tween the two theoretical models at both temperatures. 



2. Comparison with sample growth curves 

In fll"H ] some specific growth curves are also presented, and 
we shall compare these with with our computations, and those 
of Bijlsma et al. [fL7[], which appeared in preprint form after 
the initial submission of this paper. 

In Fig. [fo| we show the data from Fig. 3 of jTT|], the compu- 
tation of Fig. 1 1 of [|l7|], and our own computations. This is 
for the MIT sodium trap, with the simulation parameters taken 
from Ref. ||] of N t = 60 X 10 6 , T. t = 876nK and 77 = 2.5. 
We find this results in a condensate of 6.97 x 10 6 atoms at 
a temperature of 604nK, and a final condensate fraction of 
21.8% after half a second, which agrees wi th ou r predictions 
from the solution of the equations in Sect. V A at t = 00 to 
within 0.2%. 

We can see that there is little difference in the results of 
the two computations for this case, the main discrepancy be- 
ing that the initiation time for our simulation is a little longer 
than that of Bijlsma et al. This is likely to be due to the fact 
that their calculation starts with the condensate already occu- 
pied with ?io = 5 x 10 4 atoms, whereas we begin with the 
equilibrium number at this temperature given by the Bose dis- 
tribution of uq = 208 atoms. This difference could be brought 
about by the use of a slightly different density of states, which 
is also the likely cause of the difference in the final condensate 



number, of approximately 3 x 10 atoms. 

The agreement with the experimental growth curve data is 
very good for both computations. The simpler model of [[|] 
and QKVI cannot reproduce the results at this temperature, as 
is shown by the lower grey curve in Fig. |l0| This is as we 
expect — the final condensate fraction is far greater than 10% 
and in this case the "overshoot" of /i e s is significant. 

Given these initial conditions, this is the only case in which 
we have found that the "speedup" given by the full quantum 
Boltzmann theory may yield a significant improvement of the 
fit to the experimental data. 

We would like to emphasise, however, that the parameters 
used for this simulation do not come from Ref. jlip. The MIT 
paper does not provide any details of the size of the thermal 
cloud, or the temperature at which this curve was measured, 
and as such, a set of unique initial and final parameters of the 
experiment cannot be determined. We have simply taken these 
parameters from the calculation of Ref. [ 17 1. 



ing Fig. 5 of Ref. JIT] ] shows that most condensates of 7 x 10 6 
atoms or more were formed at temperatures above 800nK. We 
have therefore performed a second calculation using the sim- 
ple model with a final temperature of 830nK, and this result is 
shown as the upper grey curve in Fig. [l(]. As can be seen, this 
also fits the experimental data extremely well. The conden- 
sate fraction at this higher temperature is 10.2%, meaning that 
these parameters are very similar to the situation considered 
in Fig. |[ which was originally found to be a good match to 
experimental data in Ref. [^| . We note that the solution to the 
simple model at this higher temperature is also in good agree- 
ment with our more detailed calculation for these parameters. 
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FIG. 10. A comparision between the results of Fig. 11 of Jl7| ] 
and our own calculations with the initial conditions N t = 60 x 10 6 
atoms, Ti = 876nK, 77 = 2.5. Our data is shown as the solid line 
and the results from Bijlsma et al. are the dashed line. The results of 
the simple model of condensate growth with a final temperature of 
Tf = 604nK matching these initial conditions is lower grey curve. 
The upper grey curve is also for the simple model, but with what 
we feel is a more realistic final temperature of Tf — 830nK. The 
experimental data points are the solid dots. 

The situation is very different, however, if we compare with 
the result of Fig. 4 of [fit], in which the final condensate num- 
ber was 1.2 x 10 6 atoms. In this case, the data of Fig. 4(b) 
of Ref. [ pi] ] can be used to extract all the relevant experimen- 
tal parameters. This graph shows an experimentally measured 
reduction in the thermal cloud number from about 40 x 10 6 
atoms to about 15 x 10 6 over the duration of the experiment. 
Including the final condensate population gives a total num- 
ber of atoms in the system of approximately 16.2 x 10 6 , or a 
loss of about 60% of the atoms. With the three pieces of data 
taken from the MIT graphs (initial thermal cloud number, fi- 
nal thermal cloud number and final condensate number), we 
can estimate all the relevant parameters using the equations 
of Sect. V A|, and these are shown in the second column of 
Table llV 



In fact, it seems likely to us that the final temperature for the 
experimental curve shown in Fig. [l(] should be higher. Study- 



While the parameters we present here are consistent with 
the static experimental data, the growth curve corresponding 
to these parameters (shown in Fig. [TI|) certainly does not fit 
the dynamical data. We find that to remove such a large pro- 
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portion of atoms, yet still obtain a relatively small condensate, 
the initial system must be a long way from the transition tem- 
perature, with /ijnit = —IYITvJj. 

This means that condensate growth does not occur until the 
relaxation of the thermal cloud is almost complete, resulting 
in a very long initiation time, Also, when the growth does be- 
gin, the rate is significantly slower than was experimentally 
observed. This is the region in which the experimental and 
theoretical discrepancies lie. 





Parameters 


Parameters 




extracted from 


which give an 


Quantity 


experiment 


apparent fit 


iVi (10 B ) 


40.0 


40.0 


Atoms lost 


60% 


94% 


Condensate fraction 


7.2% 


51% 


Ti (nK) 


945.5 


765 


Tf (nK) 


530 


211 


V 


2.19 


0.60 



TABLE IV. Comparison of the static parameters of the Bose gas 
that match Fig. 4 of Ref. Q 

The comparison of results is presented in Fig. O. As well 
as the computation based on the extracted parameters, we also 
present two "apparent fits", one based on our calculations and 
another based on a calculation of jl7[], and here we find the 
results of the two different formulations are almost identi- 
cal. The difference appears to be due to the initial condensate 
number — our calculations begin with 295 atoms, whereas Bi- 
jlsma et al. begin with 10 4 atoms. The initial parameters cho- 
sen in [|l7| for this simulation are a system of Ni = 40 x 10 6 
atoms at at temperature of Ti = 765nK, and the energy distri- 
bution is truncated at 7] = 0.6 — an extremely severe cut. 



CO 

O 
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FIG. 11. A comparision between the data of Fig. 4 of Jll| ] (large 
solid dots) and our own calculations. The solid curve shows the 
growth curve for the static parameters that we have extracted from 
the experimental data: Ni = 40 x 10 6 , Ti = 945. 5nK, r\ = 2.19. An 
apparent fit can also be obtained — the parameters for the grey curve 
(our results) and dashed curve (Bijlsma et al.) are Ni = 40 x 10 6 
atoms, Ti — 765nK, r\ = 0.6. However, as noted in the text, these 
parameters are not experimentally acceptable. 



However, while the fit to the experimental data looks very 
good, the initial parameters for these calculations are not con- 
sistent with the experiment. An inspection of the final state of 
the gas explains the situation. The final temperature accord- 
ing to these computations is Tf = 211nK, and the condensate 
fraction is 51%. Looking at the data of JTTt], we find no re- 
ported temperatures to be lower than 500nK, and the largest 
condensate fraction reported to be 30% (although our analysis 
of their data from Fig. 5 gave a maximum of 17%). The evap- 
orative cooling of these particular simulations would have to 
remove 94% of the atoms in the trap, and we believe it is very 
unlikely that this matches any of the experimental situations. 



3. Speedup of condensate growth compared to simplified theory 

We have shown that a significant speedup of the condensate 
growth can occur at higher condensate fractions, but this can- 
not explain this particular discrepancy with the experimental 
results, for all of the measured values of temperature and con- 
densate fraction for which growth rates are presented in Fig. 5 
of 101. 

The only situation in which this speedup might possibly be 
relevant to experiment is the single growth curve correspond- 
ing to Fig. [jJJ. However, as we have noted, the initial con- 
ditions for this figure are quite speculative, and in fact also 
appear to be unrealistic. 

The actual speedup observed in our computations is of the 
order of magnitude of that achievable with a different con- 
densate fraction, and it is conceivable that the problem could 
be experimental rather than theoretical — a systematic error in 
the methodology of extracting the condensate number from 
the observed data could possibly cause the effect. For a real- 
istic comparison to be made between experiment and theory, 
sufficient data should be taken to verify positively all the rele- 
vant parameters which have an influence on the results. Thus, 
one should measure initial and final temperatures, the final 
condensate number, the number of atoms in the vapor initially 
and finally, and the size of the "cut". It should be noted in par- 
ticular, in the one case where all of this data is available — that 
presented in Fig. [n] — good agreement is not found. 

The previous paper in this series, QKVI, considered in de- 
tail a semi-classical method of fitting theoretical spatial dis- 
tributions to the two dimensional data extracted by phase- 
contrast imaging of the system during condensate growth. 
This method shows that significantly different condensate 
numbers and temperatures are consistent with the MIT data 
and methodology 1 1 1 ] . This seems to us to be a more likely 



origin of the discrepancy between theory and experiment at 
low temperatures with a small condensate number. 



H. Outlook 

It does remain conceivable, however, that approximations 
made in this formulation of Quantum Kinetic theory are not 
appropriate to the experimental regime where the discrepancy 
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remains. In this section we summarise the possible further 
extensions. 

The first is the ergodic approximation, that all levels of a 
similar energy are assumed to be equally occupied. From the 
results of QKII it would seem than any non-ergodicity in the 
initial distribution would be damped on the time-scale of the 
growth — therefore the effect of this could be significant if the 
initial distribution is far from ergodic. It is difficult to know 
what the exact initial distribution of the system is without per- 
forming a three-dimensional detailed calculation of the evap- 
orative cooling, which would require massive computational 
resources. There is also the fact that we have used the simpli- 
fied form (|l8|), derived in analogy with the work of Holland et 
al. [ pj[ ] on the ergodic approximation. 

The second important approximation is that the lower-lying 
states of the gas are reasonably well described by the single 
particle excitation spectrum, and thus using a density of states 
description in calculating the collision rates of these levels. 
The justification of this is that these states are not expected 
to be important in determining the growth of the condensate, 
and in QKVI it was shown that varying these rates by orders 
of magnitude had little effect on the growth curves. 

A third approximation made is that the growth of the con- 
densate level is adiabatic, and its shape remains well described 
by the Thomas-Fermi wave function. This may not be the 
case, and indeed some collective motion during growth was 
observed in [|PT|]. We feel this may become important for suf- 
ficiently large truncations of the thermal cloud, in experiments 
that could be considered a temperature "quench". Removing 
this assumption would require introducing a full description 
of the lower lying quasiparticle levels, and a time dependent 
Gross-Pitaevskii equation for the shape of the condensate. 

The final approximation is that fluctuations of the occupa- 
tion of the quantum levels are ignored. 

The agreement between the theory and the single experi- 
ment performed so far is generally good, and there is only one 
regime in which there is significant discrepancy. The removal 
of these approximations requires a large amount of work, and 
we feel this is not justified until new experimental data on 
condensate growth becomes available. 



VI. CONCLUSION 

We have extended the earlier models of condensate growth 
and QKVI to include the full time-dependence of both 
the condensate and the thermal cloud. We have compared the 
results of calculations using the full model with the simple 
model, and determined that for bosonic stimulation type ex- 
periments resulting in a condensate fraction of the order of 
10%, the model of [||] and QKVI is quite sufficient. 

However, for larger condensate fractions the depletion of 
the thermal cloud becomes important. We have introduced the 
concept of the effective chemical potential /i e g for the thermal 
cloud as it relaxes, and observed it to overshoot its final equi- 
librium value in these situations, resulting in a much higher 
growth rate than the simple model would predict. Thus we 



have identified a mechanism for a possible speedup that may 
contribute to eliminating the discrepancy with experiment. 

We have also found that the results of these calculations can 
be qualitatively explained using the effective chemical poten- 
tial of the thermal cloud, fag, and the simple growth equation 
(j|). In particular, the rate of condensate growth for the same 
size condensate can be remarkably similar over a wide range 
of temperatures; in contrast, the rate of growth is highly sen- 
sitive to the final condensate number at a fixed temperature. 

This model we have used in this paper eliminates all the ma- 
jor approximations in the calculation of condensate growth, 
apart from the ergodic assumption, whose removal would re- 
quire massive computational resources. In the absence of ex- 
perimental data sufficiently comprehensive to make possible a 
full comparison between experiment and theory, this does not 
at present seem j ustified. 

In Sect. V G 2 we have compared the results of our simula- 
tions to those of Bijlsma et al. ]17|], and found that our formu- 
lations are quantitatively very similar, giving growth curves in 
very good agreement with each other. The two treatments are 
based on similar, but not identical methodologies, and have 
been independently computed. Thus the disagreement with 
experiment must be taken seriously. 
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